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Abstract 

A search for double beta decays of the palladium isotopes 110 Pd and 102 Pd into excited states of their 
daughters was performed and first half-life limits for the 2^/3/3 and Oj^/3/3 decays into first excited + and 2 + 
states of 5.89 x 10 19 yr and 4.40 x 10 19 yr (95 % CL) for the 110 Pd decay were obtained. The half-life limits 
for the corresponding double electron capture transition of 102 Pd are 7.64 x 10 18 yr and 2.68 x 10 18 yr (95 % 
CL) respectively. These are the first measurements for 102 Pd. 
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1. Introduction 

During the past 20 years vast progress has been made in unveiling the properties of neutrinos. For 
decades neutrinos were thought to be massless, which no longer holds true: flavour oscillations found in the 
leptonic sector, studying neutrinos coming from the sun [HHf, the atmosphere high energy accelerators 
beams 0, Q and nuclear power plants Q , are explained by neutrino oscillations requiring a non zero neutrino 
mass. However, no absolute mass scale can be fixed with experiments studying the oscillatory behaviour. 
To achieve this, one has to investigate weak decays, such as beta decays or neutrinoless double beta decays 

(Z, A) -+{Z + 2, A) + 2e~ (0^-decay) . (1) 

The latter violates total lepton number by two units and thus is not allowed in the Standard Model. 
The 0^/3/3-decay is the gold plated process to distinguish whether neutrinos are Majorana or Dirac particles. 
Furthermore, a match of helicities of the intermediate neutrino states is necessary which is done in the 
easiest way by introducing a neutrino mass. This mass is linked with the experimentally observable half-life 
via 



(l?/ a ) * = G^(Q, Z) \M°? T - M°*f 



(2) 



where (m Ve ) is the effective Majorana neutrino mass, given by (m Vt .) — |X^^ei TO «| ana - is the 
corresponding element in the leptonic PMNS mixing matrix, G 0v {Q, Z) is a phase space factor and M*q T — 
Mp u describes the nuclear transition matrix element. The experimental signature is the emission of two 
electrons with a sum energy corresponding to the Q-value of the nuclear transition. Apotential evidence 
has been claimed in the 0j//3/3-decay of 76 Ge with T°J 2 = 2.23±g;^ x 10 25 yr at 90% CL Q. In addition, the 
SM process of neutrino accompanied double beta decay, 

(Z, A) -> (Z + 2, A) + 2er + 2v e (2^-decay) (3) 
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can be investigated, which is expected with half- lives arround 10 20 yr. For recent reviews see [8|. 
Additional information is provided by the alternative process of positron decay in combination with 
electron capture (EC). Three different decay modes can be considered: 

(Z,A)^(Z-2,A) + 2e+ + (2v e ) (/J+/3+ ) (4) 

e~ + {Z,A) -^{Z -2,A) + e+ + {2u e ) (/3+ -/EC) (5) 
2e~ + (Z,A) -+{Z -2,A) + {2v e ) (EC/EC) (6) 

Decay modes containing a positron have a reduced Q-value as each generated positron accounts for a 
reduction of 2 m e c 2 . Thus, the full energy is only available in the EC/EC mode and makes it the most 
probable one. However, it is also the most difficult to detect, only producing X-rays instead of 511 keV 
gammas. Furthermore, it has been shown that /3 + /EC transitions have an enhanced sensitivity to right- 
handed weak currents (V+A interactions) [9j and thus would help to disentangle the physics mechanism of 
O^/3/3-decay. In the last years, also neutrinoless EC/EC modes have been discussed with renewed interest, 
because of a potential resonance enhancement up to a factor of 10 6 in the decay if the initial and final 
excited state are degenerate [lfj . Recently, a series of isotopes with extreme low Q-value was explored for 
enhancement in the ECEC mode and with 152 Gd a very promising candidate was found (Tlj . 

Another branch of search is linked to excited state transitions. The signal in this case is extended by 
looking at the corresponding de-excitation gammas. However, in the approach of a passive sample on a 
Ge-detector it will not allow to distinguish between the 2z//3/3-decay and O^/3/3-decay mode. Thus, the 
deduced half-live limits are valid for both. The investigation of 2^/3/3-decay modes into excited states will 
add information on nuclear structure, valuable for matrix element calculations. Furthermore, a potential 
observation of Oz^/3/3-decay into an excited 2 + -state would likely point to other contributions besides neutrino 
masses. The searches described in this paper are based on the search for excited state transitions. 

An element getting little attention in the past is palladium with the isotopes of interest 110 Pd and 102 Pd. 
Among the eleven Oi//3/3-decay candidates with a Q-value larger than 2MeV, 110 Pd has several advantages: 
it has the second-highest natural abundance (11.72%) and, in addition, it is an excellent candidate to 
probe the single-state dominance hypothesis for 2i//3/3-decay, i.e. that only the lowest lying intermediate 
l + -state will contribute to the nuclear transition matrix element describing its 2i^/3/3-decay. Only one rather 
weak experimental limit in the order of 10 17 yr exists for Oi//3/3-decay ground state transitions in 110 Pd [HI ]. 
Theoretical predictions for the 2^/3/3-decay ground state mode are in the range of 0.12 - 29.96 x 10 20 yr 
I3I 14, Til 16, 17. fl8l Til 21 1 . Theoretical predictions for the excited state transitions are 4.4 x 10 25 yr 



8.37 x 10 25 yr [lfj, 1.5 x 10 25 yr [20j and 0.62 - 1.3 x 10 25 yr [2l| for the 2+ state and 2.4 x 10 26 yr 



and 4.2 - 9.1 x 10 23 yr [2l| for the 0^ state. 
The second isotope 102 Pd has a Q-value of 1172 keV, a natural abundance of 1.02% and is able to decay 
via EC/EC and /3 + /EC. These decay modes have never been studied for 102 Pd experimentally and no 
theoretical predictions exist. The level schemes of both isotopes are shown in Fig. [1] 



2. Experimental Setup 

The measurement was performed in the Felsenkeller Underground Laboratory in Dresden with a shielding 
depth of 120 mwe. A sample of 802.35 g of Pd was used, which was purified before by C. HAFNER GmbH 
+ Co. KG. It was placed in a standard Marinelli baker (D6) of 70 cm diameter and 21cm height. The 
sample was positioned on a HPGe detector with an efficiency of 90 % routinely used for 7-spectroscopic 
measurements. The detector has an 1 mm thick Al-window towards the sample. It is surrounded by a 
5 cm copper shielding embedded in another shielding of 15 cm of clean lead. The inner 5 cm of the lead 
shielding have a low contamination of 210 Pb of only 2.7 ± 0.6Bq/kg while the outer 10 cm have an activity 
of 33±0.4Bq/kg. 

The detector is placed inside a special building which acts as a Faraday cage and a Rn shield; additionally, 
the setup is flushed with nitrogen in order to reduce Rn contamination. More details can be found in [22I [23| 
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Figure 1: Level schemes of 110 Pd (left) and 102 Pd (right) decays. 

and a schematic drawing is shown in Fig. [5J The Pd was stored underground for more than one year prior 
to the measurement except for 18 d of purification. The data were collected with a 8192 channel MCA from 
Ortec and were converted into the ROOT format for analysis. 
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Figure 2: Schematic drawing of the used setup. Composed from pictures in [23 



An extensive calibration was performed using 17 7-lines from 
energy calibration curve of 



different nuclides resulting in a linear 



£[kcV] = 0.342746 ■ 



keV 



x channel - 4.337734 keV. 



(7) 



channel 

The measurement range of the spectrum goes up to 2.8 MeV. The resolution was calibrated using the 
calibration lines and fitted with a second order polynomial. The actual values at the energies of the lines of 
interest will be discussed in the corresponding analysis section. 

Despite purification, the measured spectrum is dominated by intrinsic contaminations of the Pd. Clear 
7-lines from the 238 U and 232 Th decay chains as well as 40 K are visible. However, former Americium 
contributions have been removed completely by the purification. The actual background spectrum of the 
detector system itself without any sample is at least an order of magnitude smaller in the regions of interest 
and can be neglected, hence the spectrum is completely dominated by the Pd sample contaminations. 
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3. Analysis and Results 



A total of 16.2 d of data were accumulated resulting in 13.00 kg • d of exposure. In the following the two 
isotopes of interest are discussed separately. For the analysis, the natural abundances of the latest IUPAC 
evaluation have been used which are 11.27% ( 110 Pd ) and 1.02% ( 102 Pd) respectively [24j. As the searches 
are purely based on gamma detection, the obtained results apply for both, Oz//3/3-decay and 2i//3/3-decay 
modes. The decays into the first excited 0^~-state de-excite via an intermediate state. Thus, there 
will be an angular correlation among the gammas, with an observational probability W(0) that the second 
gamma is emitted with an angle 9 with respect to the first one given by 



It can be seen that the probability of both gammas bccing emitted in the same direction is larger than 
an uncorrelated emission. However, the detection efficiency for a single gamma is small which results in 
a low probability to observe a summation peak and thus the searches are based on the individual gamma 
energies only. 

The efficiency for full energy detection was determined and cross checked in several ways. The most 
important one was replacing the actual used volume by a SiC>2 sample of exactly the same geometry. The 
intrinsic contaminations of the natural decay chains of 238 U and 232 Th as well as 40 K produced various 
7-lines and acted as an extended calibration source. The 17 7-lines used for the energy calibration were also 
used for the efficiency determination in the region from 238.6 keV (from 212 Pb decay) up to 2614.3 keV (from 
208 T1 decay). The efficiency in the region above 200 keV can be well fitted by two exponential functions. 
It varies between 5.77% at 468 keV and 3.89% at 815 keV, being the lowest and highest energy lines of 
interest for the search described in this paper. To account for the difference in self-absorption of Pd and 
SiC>2, extensive Monte Carlo simulations were performed using the AMOS code (25|. The amount of Monte 
Carlo was chosen in a way that the statistical error for the full energy peak for both, Pd and Si02, was 
less than 0.1%. The simulations agree within an error of less than 15% with the measurements and tend 
to be slightly higher. This can easily be explained by small geometric differences in the simulation and the 
experiment. However, independent from that is the ratio of both self-absorption simulations for Pd and 
Si02- Hence, the ratio was used to scale the experimentally well determined efficiency curve of Si02 to the 
one of Pd (Fig. [3]). A validation of the procedure was performed within the vicinity of the lines of interest 
by using prominent background lines apparent in the spectrum, namely 295.21 keV and 351.92 keV (from 



Pb), 233.63 keV (from 212 Pb), 583.19keV (from 208 T1) and 609.32keV (from 214 Bi), which results in good 



The analysis is based on the extraction of upper count limits with the Feldman Cousins method described 
in (2(|. A simple constant was chosen as a background model after concluding that the area around the 
ROI is sufficiently flat. The background was determined with a likelihood fit in a selected region around the 
peak position taken to be ±30keV excluding the peak range which is considered to be ±5keV. 

Likelihood fits with the Gaussian peak shapes and the background constants resulted in statistical 
downward fluctuations for all four peaks. This is a clear sign that no significant signal is observed (see Fig. U 
and[S]). In this case, an upper limit for the count rate can be calculated using the background only hypothesis 
and the maximal statistical fluctuation of a Gaussian distributed background for a certain confidence level. 
This is commonly referred to as sensitivity. However, the Feldman Cousins approach also considers the 
observed downward fluctuations and results in more appropriate results for low count rates. 

In order to obtain a numerical value for the upper limit, all bins within the FWHM of an expected peak 
are combined into a single analysis bin. The resulting fraction of the peak that is covered by the analysis bin 
is dependent on the binning of the data but always larger than 76 %. The background expectation and the 
measured count rate are used to evaluate an upper count limit for this bin with the ROOT implementation 
of the Feldman and Cousins method. The result is then scaled to the full peak area. 

From the background point of view, exactly the same 7-lines in the performed search could be produced 
from the beta decays of the intermediate nuclide of the investigated double beta system, which will be 




(8) 



agreement. 
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Figure 3: Measured efficiencies and the fitting function. 



discussed in more detail in the next section. However, their contribution can be rejected by the non- 
observation of other, more prominent 7-lines at different locations in the spectrum. The only prominent 
background line to be expected within ±5 keV of any of the four lines under investigation is from 137 Cs at 
661.66keV potentially influencing the 2^-limit in the 110 Pd system. However, no indication of this line is 
observed. 

3.1. The 110 Pd system 

Two lines were investigated for 110 Pd at energies 657.76 keV (corresponding to the 2f transition) and 
815.35 keV (additionally emitted in the 0^ decay) respectively. The corresponding energy resolutions at 
these energies are 1.51 and 1.61 keV (FWHM). Potential 7-lines mimicking the signal would come from 
110 Ag and 110m Ag decays. 110 Ag has a half-life of 24.6 s only and thus has to be produced in-situ. With 
the given shielding this can be excluded. Potentially more dangerous is the long-living 110m Ag (half-live of 
249.79 d). This isotope has two prominent lines at 1384.3 keV and 1505.04 keV. They are not visible in the 
spectrum and thus this contribution can be excluded for this search. No lines were visible at both peak 
positions of interest and thus an upper limit (95 % CL) of signal events of 10.53 and 7.34 could be extracted 
for 657.76 keV and 815.35 keV respectively. 

Using the known Pd mass and efficiencies, this can be converted into lower half-live limits of 



T (0v+2 V ) 110pd ^ no Cd ( + 8 i5.3 fey) > 5.89 x 10 19 yr (95%CL) (9) 
T (0v+2 V ) 110pd ^ iio Cd ( 2 + 657.8 keV) > 4.40 x 10 19 yr (95%CL) . (10) 

These are the first experimental limits for excited state transitions in the 110 Pd system. 
3.2. The 102 Pd system 

Two lines were investigated for 102 Pd at energies 468.59 keV (only emitted in the Of decay) and 475.05 keV 
(corresponding to the 2f transition) respectively. The energy resolution at these energies is 1.39 keV 
(FWHM) for both lines. Potential 7-lines mimicking the signal would come from 102 Rh and 102m Rh decays. 
102 Rh with a half-life of 207 d has no reasonable line to check. The strongest one is the 475.05 keV line. As 
there is no signal in this region it can be concluded that it does not contribute to the 468.59 keV region. 
On the other hand 102m Rh (half-live of 2.9 yr) has multiple lines to explore, the most restricting ones are 
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a line at 631.28 keV with 56% emission probability and at 697.49 keV with 44%. Both of them are not 
observed in the spectrum and thus can exclude such a contribution. No lines were visible at both peak 
positions of interest and thus an upper limit (95 % CL) of signal events of 17.64 and 6.24 could be extracted 
for 475.05 keV and 468.59 keV respectively. The obtained half-live limits are 

r (OH-2«0 102pd ^ i02 Ru(0 + 468 6keV) > 7 . 6 4 x 10 is yr ( 95 %CL) (11) 
T (°v+2v) io 2pd i02 Ru ( 2 +, 475.1 keV) > 2.68 x 10 18 yr (95 %CL) . (12) 
These are the first experimental limits on 102 Pd double beta decays. 



4. Summary 

Double beta decay transitions into excited states for the two Pd-isotopes 102 Pd and 110 Pd have been 
investigated for the first time. These transitions contain valuable informations about the physics mechanism 
of double beta decay and the involved nuclear physics. However, no signal into the first excited + and 2+ 
states have been observed. 
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Figure 4: Peak region of 102 Pd for the + — > 2 + transition with the illustration of the background model, the best fit 
and the peak with the upper count limit at 95 % CL. 
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Figure 5: Peak region of Pd for the + — > 2 + transition with the illustration of the background model, the best fit 
and the peak with the upper count limit at 95 % CL. 
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